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Abstract. Aims: We investigate the behaviour of dust in protoplanetary disks under the action of gas drag in the presence of a 
planet. Our goal is twofold: to determine the spatial distribution of dust depending on grain size and planet mass, and therefore 
to provide a framework for interpretation of coming observations and future studies of planetesimal growth. 
Method: We numerically model the evolution of dust in a protoplanetary disk using a two-fluid (gas + dust) Smoothed Particle 
Hydrodynamics (SPH) code, which is non-self-gravitating and locally isothermal. The code follows the three dimensional 
distribution of dust in a protoplanetary disk as it interacts with the gas via aerodynamic drag. In this work, we present the 
evolution of a minimum mass solar nebula (MMSN) disk comprising 1% dust by mass in the presence of an embedded planet. 
We run a series of simulations which vary the grain size and planetary mass to see how they affect the resulting disk structure. 
Results: We find that gap formation is much more rapid and striking in the dust layer than in the gaseous disk and that a system 
with a given stellar, disk and planetary mass will have a completely different appearance depending on the grain size. For low 
mass planets in our MMSN disk, a gap can open in the dust disk while not in the gas disk. We also note that dust accumulates 
at the external edge of the planetary gap and speculate that the presence of a planet in the disk may enhance the formation of a 
second planet by facilitating the growth of planetesimals in this high density region. 
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1. Introduction 

The effect of a planet in a gaseous di sk has been well stud- 



ied both analytically and numerically dGoldreich & Tremaine 



Lin & Papaloizou 





1993j 


Kley 





1999; Bryden et al. 




de Val-Borro et al.l 120061) . Planetary gaps result from tidal in- 



teractions between a sufficiently massive planet and the disk 
if the tidal torques induced by the planet exceed the vis- 
cous torques of the disk. If the gap is sustained, its pres- 
ence severely decreases accretion of disk m aterial by the 



plane t (.Artymowicz & Lubow..l996: iKlevlll999: Bryden et al 



1999h . w hich suggests t hat the planet formation process is self- 



limiting. iLufkin et al 



(l2004 l) suggest, however, that the ac- 
cretion rate is not so dramatically reduced in self-gravitating 
disks. The exchange of angular momentum that results from 
the tidal torques leads to planetary migration, which is slowed 
once a gap has formed. This is kno wn as Type II migration 



and proceeds on a viscous timescale dPapaloizou & Lin 1 1984 
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Takeuchi et al. 119961: I WardI 119971) . If the planet mass is not 
massive enough t o sustain a gap, it will unde r go the faste r 



Type I migration ( Goldr eich & Tremaine 119791: 



Ward 1997) 



while intermediate-mass planets only o pen a partial gap and 
suffer Type III, or runaw ay, migration (iMasset & Papaloizou 
20031: lArtymowiczll2004 . 



There are a variety of observational signature of planetary 
gaps, including infrared dips and direct scattered light and sub- 
millimetre observations of protoplanetary disks. Mid-infrared 
dips are seen in quite a few pre-main sequence disks around 
T Tauri stars including T Tau and GM Aur The mid-IR emis- 
sion comes from the inner 1-10 AU of the disk and the 10 /urn 
feature is predominantly from silicate grain emission and it 
has been suggested that such IR dips may be the results of an 



unseen planet clearing out the disk m aterial (e.g. lCalvet et al 



20021: iRice et all 120031; IVarniere et al.ll2006.) . Such dips, how 



ever, have also been modelled assuming a continuous disk with 
a more realistic dens ity distribution than a standard power-law 
dBoss & Yorke 1996). and by a continuous disk with silicate 
grain growth dDuUemond & DominikI l2005b . Sub-milUmetre 
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observations of older debris disks show brightness asymme- 
tries and a range of structures within the disk, many of which 
can be modelled by th e presence of a n unsee n planet locking 
the dust in resonances ( Ozernov et al.L 2000; Wil ner et al.l2002t 
Deller & Maddisonl2005b" H oweve r it has been noted bv lWvattI 
(l2006h and iGrigorieva et alJ (l2007h that holes and dumpiness 
of debris disks can result pu rely from coUisional cascade of 
planetesimals within the disk. lTakeuchi & Artymowica (1200 ih 
show that gas drag alone can be used to explain annulus-like 
disks with sharp edges, and sub-mm holes have als o bee n 
explained by the sublimation of icy grains djura et allll998l) . 
Direct imaging of scattered light from protostellar and debris 
disks with HST confirm the presence of g aps (Weinberaer et al. 
1999; S chneider et all 1 1999t lArdila et al l |2004 iKalas et all 
2005t , Schneider et al.ll2O06r 

Recent models have indicated that ALMA will be able to 
detect planetary gaps at sub-millimetre wav elengths. Using the 
results of 2D hydrodynamics simulations. Wolf & D'Angelol 
(120051) use a Monte Carlo radiative transfer code to predict 
the emission of a disk hosting a planetary gap carved out by 
a range of planet masses. They find that ALMA will be able 
to detect the gap in a O.OIMq disk with an Mjup planet at 
a distance of 100 pc. When using 2D hydrodynamics codes, 
one needs to add an analytic disk atmosphere in order to pro- 
vide the Monte Carlo simulation with a 3D density struc- 
Wolf& D'Angelol (E005l) me an un-flared disk, while 



ture. 



Varniere et al. ( 200 61) use a flared disk and show that the spec 
tral energy distribution and synthetic scattered light images are 
effected by the outer wall of the gap. Both groups, however, 
assume that the gas and dust are well-mixed within the disk 
and yet the gas-to-dust ratio will change substantially in pro- 
toplanetar y disks as the grains b egin to grow and settle to the 
midplane (IMaddison et al ]l2003h . 

In recent years, a growing number of studies started to 
consider the separate evolution of gas and dust in disks. 
Barriere-Fouchet et al.l (12005 , hereafter BF05) described the si- 
multaneous radial migration and vertical settling of dust grains, 
treated as a separate fluid in glob al 3D SPH si mulations, 
and their dependency on grain size. iRice et al.l (12004) showed 
that large solid bodies, treated as trace particles, tend to ac- 
cumulate in the spiral arms in massive di sks. The influence 
of the magneto-rotational instability (MRI, iBalbus & Hawlev 
199ll) ha s also been studied in E ulerian global 2D sim- 
ulatio ns (jproman g & Nelson 2005) and local 3D simula- 
tions (jpromang & Papaloizou ,2006: iJohansen & Klahi 2005 



Johansen et al.l2006l) . A fraction of the sohd bodies, treated 



as a separate fluid for small sizes or as particles for large sizes, 
were found to be trapped in density maxima generated by the 
MHD turbulence. 

As well as the gas-to-dust ratio varying throughout the disk, 
we expect that the effects of planetary gaps will be stronger 
in the dust phase than in the gas phase, which will further af- 
fect observations. Any gap created by a planet will result in a 
pressure minimum in the disk, and unless the dust and gas are 
extremely well coupled and the motion is governed by the gas 
(which will be true for sub-micron sized grains), the gas will 
respond to the pressure gradient while the dust will not. (The 
dust feels no viscous forces and much lower pressure forces 



than the gas and hence the tidal torques which create the gap 
will dominate.) If anything, the dust moves to regions of pres- 
sure maxima ('see iHaghighipour & Bossll2003 ) on the edge of 
the gap. Further more, it is e asier to form gaps in 2D disks than 
in 3D disks (see iKIevll 19991) and we expect the dust to settle to 
the midplane on a much shorter timescale than the gas phase. 
Both of these phenomena can occur at the same time. Thus 
predictions of observations of gap which assume a well-mixed 
disk are probably too pessimistic and it is therefore important 
to study the effects of a planet on the dust layer to help con- 
strain the characteristics of the planets we might detect via the 
gaps they create. 

While many semi-analytic and numerical studies have in- 
vestigated the formation of planetary gaps in gas disks, very 
little work had been done on the effect of a planet in a dusty 
disk. Over the last few years, however, attention has begun to 
focus of the du st phase of 2D protoplanetary disks wit h em- 
bedded planets. jPaardekooper & MeUemal (|2004 12006 !) stud- 
ied the conditions for gap formation and showed that a smaller 
planetary mass is required to open a gap in the dust disk rather 
than in the gas. They note, however, that one should be care- 
ful when interpreting results f or planets l i ghter than O.lMjup in 
two dimensional simulations. Rice et al. I (120061) found that the 
pressure gradient at the outer gap edge can act as a filter keep- 
ing large grains out of the inner disk, which may explain the 
deficit of near-IR flux observed in the SED of some T Tauri 



stars. [Alexander & Armitagd (l2007l) studied the dynamics of 
dust and gas during disk clearing and, adding an embedded 
planet, show that observations of transition disks can help dis- 
criminate between different models of disk clearing. 

In this paper we study the formation of a gap triggered 
by an embedded planet in the dust layer of a protoplanetary 
disk. We will study the effects of planet mass and grain size 
on gap formation and evolution in 3D, two-phase (gas+dust) 
protoplanetary disks. These result will be compared with the 
gas-only simulations of others who must then assume the gas 
and dust are well mixed in order to determine how the dust- 
gas coupling (or decoupling) affects potential observations of 
planetary gaps. 

2. Gap opening criteria 

The gravitational perturbation of the planet launches a spiral 
density wave at the Lindblad resonances of the disk. The re- 
sulting tidal torques lead to an exchange in angular momen- 
tum when the waves damp by viscous dissipation or shocks. 
Waves exterior (interior) to the planet carry excess (deficit) an- 
gular momentum compared to the disk and hence when the 
wave damps, the disk gains (loses) angular momentum and 
moves outwards (inwards). This results in the formation of a 
gap around the planet. 

In order to sustain the gap, there needs to be a balance be- 
tween the tidal torques, which act to clear the gap, and the vis - 
cous torques, which act to fill the gap (Lin & Papaloizou 1979*). 
The gap criterion dPapaloizou & Lia,1984>:,Bryden et al..,1999.) 
is given by 

^tidal ^ ^visc 
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> 40ffss - 



H 



(1) 



wh ere Mp and M» are the mas s of the planet and star, ass is 
the IShakura & SunvaevI (Il973h viscosity parameter, H is the 
disk scale height and is the semi-major axis of the planet. 
For inviscid disks, the competing forces are local pressure and 
gravitational tidal perturbations dLin & Papaloizoull 1993b and 
gap criterion is given by 



r tidal 



press 



^>3 ^ 
Irp 



(2) 



Our previous simulations (BF05) show that the thickness 
of the dust layer depends on the grain size because different 
sized grains fall to the midplane at different rates. The gas 
drag ensures that (for certain grain sizes which depend on the 
nebula conditions) the dust rapidly settles to the midplane and 
since the above expressions for the gap criteria depend on the 
disk scale height, this would suggest that it is easier to create 
and sustain a gap in the dust layer than in the gas. However, 
Eqs. ([T]l and (|2]i were derived without taking the drag force into 
account. Very small grains, which feel the strongest drag, are 
basically co-moving with the gas and are not likely to create 
a much sharper gap than the gas. Conversely very large grains 
(or boulders) are essentially decoupled from the gas and, feel- 
ing very little drag, have very long settling times. In this pa- 
per, we concentrate on the more interesting intermediate sized 
grains with the most efficient settling rates, for which the sit- 
uation is greatly complicated by the drag forces and therefore 
numerical simulations are necessary to fully describe their be- 
haviour The regi me we consider here is somewh at different 
to the 2D work of Paardekooper & M ellema ( 2004') who study 
dusty disks in which the dust is strongly coupled to the gas and 
responds indirectly to the planet gravity through small radial 
pressure gradients in the gas disk that leads to the formation of 
a gap. 

3. Code description 

We have developed a 3D, two-fluid (gasH-dust), locally isother- 
mal, non-self-gravitating code based on the Smoothed Particles 
Hydrodyn amics (SPH) algori thm, a Lagrangian technique de- 
scribed by iMonaghan (Il992h . The dusty-gas is approximated 
by two inter-penetrating fluids that interact via aerodynamic 
drag. We assume Epstein drag, which is valid for the range of 
dust sizes and nebula parameters used in this study. Thus the 
drag force Fd is given by 



An 



-PgS vc. 



(3) 



where pg is the gas density, s is the (spherical) grain radius, 
c is the sound speed, and v = |vd - Vg| is the velocity differ- 
ence between dust and gas, given by Vd and Vg respectively 
dStepinski & ValageasI Il996h . The dust particles are incom- 
pressible and there is no coagulation or evaporation of grains. 
All simulations in this paper are conducted with just one grain 
size at a time. In future work we shall include multiple grains 



sizes per simulation. For details of how the equations of mo- 
tion and the density of the two fluids are calculated, we refer 
the reader to BF05. 

3.1. Implementation and initial conditions 

We set up a disk of gas and dust with a total mass of Mjisk 
around a one solar mass star into which we embed a planet of 
mass Mp at a distance of rp. The system is evolved for about 
100 planetary orbits. 

The code units are chosen such that G = = rp = 1. 
The disk equation of state is isothermal and hence the temper- 
ature is constant in the vertical direction but follows a radial 
temperature profile (T oc r '). The isothermal sound speed at 
r = I is c - H/r - 0.05, where H is the disk scale height and 
r the radial distance from the star. Any heat produced in the 
disk due to tidal forces or viscous dissipation is assumed to be 
immediately radiated away. In order to compare our results to 
those of Paardekooper & Mellema (2006), we use an initially 
constant surface density profile. In the z direction particles are 
randomly distributed between +H. The intrinsic density of the 
dust grains, p^, is 1.25 g cm""*. 

All simulations contain an equal number of gas and dust 
particles and the disks contain 1% dust by mass. We start our 
simulations with a gas disk containing 50,000 particles and an 
embedded planet. This is evolved for 8 planetary orbits before 
we add the dust phase. The dust particles are initially overlaid 
on the gas particles and the system is then evolved for a further 
96 orbits, making 104 orbits in total. Given that some accretion 
occurs during the first 8 orbits, we generally end up with about 
96,000 particles in total. While this may seem to be quite low 
resolution for to day 's SPH simulations (c.f. A^gas ~300,000 in 
Rice et alj2603l and lSchafer et al.ll2004 BF05 showed that this 
level of resolution gives the same results as simulations with 
almost twice as many particles (in the case of a disk without a 
planet). It should be noted that in our two-phase code, the time 
step is severely limited by the drag timescale. While the drag 
term is solved implicitly to alleviate the cripplingly short dtdmg 
issue, this requires an extra N/2 operations per time step (and 
an additional SPH particle neighbour search) and thus particle 
resolution suffers as a consequence. To ensure the validity of 
our low resolution results, we ran a series of high resolution 
simulations with A^~400,000 which are presented in Sect. 14.11 

The planet is implemented as a point mass particle which 
moves under the gravitational influence of the central star, but 
does not feel pressure, viscous or drag forces. The SPH parti- 
cles, on the other hand, feel the gravitational potential of the 
planet. The planetary orbit can either be circular, elliptic or 
even hyperbolic and is determined by the choice of the initial 
velocity. In this work we only study circular orbits. Particles 
are accreted by the planet if they cross a critical radius given 
by 



■ O.IR 



Hill 



1/3 



(4) 



We take no account of any mass or angular momentum transfer 
corresponding to this infall and therefore do not simulate mi- 
gration of the planet. The central star also accretes particles that 
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Table 1. Parameters of our simulation suite. 



MMSN models 



Star mass 


1 Mo 


Inner radius 


2 AU 


Disk mass 


2.9 10"-'' Mo 


Outer radius 


20.8 AU 


s (m) 


Mp (Mjup) 


a, /? 


Drag 


1.0 


1.0 


0.1,0 


yes 


0.1 


1.0 


0.1, 


yes 


0.01 


1.0 


0.1,0 


yes 


1.0 


0.5 


0.1,0 


yes 


1.0 


0.2 


0.1, 


yes 


1.0 


0.1 


0.1,0 


yes 


1.0 


0.05 


0.1,0 


yes 


1.0 


1.0 


0.1,0.5 


yes 


1.0 


1.0 


1, 2 


yes 


1.0 


1.0 


0.1, 


no 


1.0 


1.0 


0.1,0 


yes 



get closer than = 0.4 code units, which sets the inner edge of 
the disk. This rather large value was chosen t o co mpare with the 
results of IPaardekooper & Mellemd (|2004 E0O6,) . (See, how- 
ever. Sect. 15. II for a discussion about the location of the inner 
boundary and its effect on the gap.) Particles are also removed 
from the simulation if they escape to large radii, which is set to 
Tesc = 4 code units. 



3.2. Simulation Suite 

In this work we present the results of a small, low mass disk 
close to the minimum mass solar nebula (MMSN). The disk 
mass is set to 2.9 10"^ Mq. Given the code scahng, Rmn = 
0.36 AU, the stellar accretion radius is ri,, - 2 AU, and 
resc - 20.8 AU, with the latter two setting the radial extent of 
the disk. Our standard model has a 1 Jupiter mass (Mjup) planet 
on a circular orbit of radius 5.2 AU in a disk containing grains 
1 m in size. Such a large grain size was chosen for our standard 
model because for these nebula conditions this is the grain size 
that settles most efficiently and thus the effect of the planet is 
most striking in a dust layer of 1 m sized "boulders". Further, 
boulders of this size should radially migrate most rapidly and 
here we will show that this migration can be stopped by plane- 
tary gaps. While the 1 m grains are observationally uninterest- 
ing, a large population of boulders would result in a population 
of smaller grains by collisions which are detectable. 

We start by running our standard MMSN model and com- 
pare the evolution of the planetary gap in the gas and dust 
phases. We then run a series of experiments to study the ef- 
fect of grain size in the dust disk, with s - I cm, 10 cm and 
1 m. This is followed by a series of experiments that study the 
effects of planetary mass on gap formation and evolution, with 
Mp = 0.05,0.1,0.2,0.5 and 1.0 Mjup. The parameters of our 
simulation suite are summarised in Table [T] 



3.3. Dust Stopping time 

Since the drag force acting on the grains depends on the nebula 
gas density and relative velocity as well as the grain size, it 
is useful to know the stopping time for different grain sizes in 



Table 2. Percentage of particles lost due to accretion or ejection 
for various SPH viscosity parameters a and /3. 



viscosity parameters 


% of particles 


a p 


lost 


0.1 


16.1 % 


0.1 0.5 


16.6 % 


I 2 


32.5 % 



10' 



>. 
■55 

c 

"D 



10" 



0) 

o 
CO 



10"' 



10" 



a=0.1, p=0 
a=0.1, p=0-5 
a=1 , p=2 
gas 




10 

r(AU) 



15 



20 



Fig. 1. Comparison of the dust disk surface density after 104 
planet orbits for a range of SPH viscosity parameters, a and yS, 
for our standard model with a 1 Mjup planet at 5.2 AU. The gas 
surface density is scaled by a factor of 0.01 for direct compari- 
son with the dust phase. The vertical lines are, from left to right, 
the 1 :2 and 2:3 internal resonances, the 1 : 1 planetary orbit, and 
the 3:2 and 2:1 external resonances respectively. 

order to compare with different n ebula models. The stopping 
time dStepinski & Valageas Il996h in the Epstein drag regime 
(assuming solid spherical grains) is given by 



PiS 



(5) 



where pd is the intrinsic dust density. This is made non- 
dimentional via fs - Tf^/Qj^sp, and thus the non-dimensional 
stopping time, Ts, can be written as 



(6) 



where Qkbp = I'Kep/'" is the Keplerian angular velocity. We find 
that at the location of the planet for our standard model (see 
Sect. 14.11 ) the stopping time is given by - 0.067icm where 
Scm is the grain size in cm. Thus - 0.067,0.67 and 6.7 for 
1 cm, 10 c m and 1 m grai ns respective ly. Thi s is in good agree- 
ment with Paardekooper & Mellernal ( 20061) (noting that they 
have an additional facto r of Vtt/S in their ex pression of the 
stopping time, following Takeuchi & Linll2002 ). 



3.4. Viscosity 

Since protoplanetary disks are viscous accretion disks, 
we can apply the Shakura-Sunyaev a-disk prescription 
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( Shakura & Sunvaevlll973 ) in which the kinematic viscosity is 
given by 



V - asscH . 



(7) 



In the thin disk approximation, Hjr - c/vKep and the viscos- 
ity can be rewritten as v = a^^c^ jO., and orss is close to 0.01 
(whic h provides a good fit to observations of protostellar disks 



■ see 



Hartmann et al. 1998' King et al 2007) . 



We use, as in BF05, the standard SPH artificia l viscosity 
term (|Monaghan & GingolJl 9831: iMonaghanll 19891 ) that is in- 
cluded in the SPH momentum equation and acts as a viscous 
"pressure". Using the SPH particle subscript notation / and j, 
the viscous term of the momentum equation is given by 



dt 



(8) 



where v is the particle velocity, m the particle mass, W the 
smoothing kernel, and 11 the ar tificial viscosity term. The most 



common form of 11, j is due to iMonaghan (.1989i) and is given 
by 



n,7 = - 



acij^iij +/3fil 



PU 



(9) 



where fiij = /iv,j ■ rij/{rfj + rf-) if \ij ■ rij < and is zero oth- 
erwise, pij and Cij are the average density and average sound 
speed of interacting particles / and j, and rj^ = 0.0 l/i^ is a soft- 
ening parameter that prevents fi/j becoming singular, and h is 
the SPH smoothing length. We use the cubic spline kernel of 
Monaghan & Lattanzio (1985) given by 



W(r,/,)=^. 



1 - lir/h? ■ 
0, 



ifO<(r/h)< 1; 

if 1 < (r/h) < 2; (10) 

otherwise. 



where D is the number of dimensions, and cr is a normalisation 
constant given by 2/3, IQ/ln and l/;r in one, two and three 
dimensions respectively. 

There are two free SPH viscosity parameters, a and 
which control the strength of the viscosity. (Note that the 
Shakura-Sunyaev ass and the SPH a are not the same.) The 
viscosity associated with the linear a term of n,j was in- 
troduced to_2emove_subsoni£veloc^ oscillations that follow 
shocks ('Monagh an & Gingold 1983 ). while the nonlinear /3 
term damps hi gh Mach number s hocks and prevents particle in- 
terpenetration ( Monaghanll989h . While it is reasonable to set /? 
to zero since we do not expect strong shocks in our protoplan- 
etary disk simulations (c.f. BryjleiLet al. 1999), some artificial 
bulk viscosity via the a term is needed to prevent random ve- 
locities from growing and producing a thick disk. Substituting 
Eqs. (|9]l and ( fTOt into Eq. dSJ, it can be shown that a is related 
to the viscosity via 



Vsph = d-ach , 



where o" = 1/8 for a 2D cubic spline kernel or 2/1 5 for a 3D 
cubic spline kernel (lPongracidll988t iMurrav Il99d) . Equating 



Eqs. (|7]l and (fTTl i. we can relate ass and a for our 3D simula- 
tions via 



id-hn\ 

ass - a \ I • 



(12) 



In order to keep the corresponding Shakura-Sunyaev viscos 
ity parameter ass close to .01 (as indicated by observations 



of pro toplanetary disks - see lHartmann et al.lll998t iKing et al 



20071), we set the a = 0.1. A problem with using SPH to model 
planetary gaps is that the numerical viscosity cannot be reduced 
to arbitrarily small values, which means that (gas) particles will 
diffuse into the gap. Many nu merical simulation of p lanetary 
gaps use grid-based codes (e.g. lD'Angelo et al.ll2002l) in order 
to keep ass as low as possible. Some grid codes can get ass 
to be as low as 10 However, since observations of proto- 
planetary disks indicate ass ~ 0.01, SPH can do an adequate 
job and indeed o ne of the first sim ulations of an embedded 



lanet in a disk bv lSekiva et all (1 19871) used SPH. ISchafer et al 



(•2004) present a new treatment of the viscosity in SPH codes 
which includes an additional artificial bulk viscosity term in the 
Navier-Stokes equation and demonstrate that their results com- 
pare well with grid-based codes. This is however very compu- 
tationally expensive and has not been implemented in our code. 

To check the validity of the SPH viscosity parameters, we 
ran a series of tests to investigate the effect of changing a and yS. 
Our standard model has a = 0.1 and jS = 0, and we also ran test 
simulations with (a,y6)=(0. 1,0.5) and (1,2). As expected, we 
find the (a,y6)=(l,2) case to be much more dissipative than the 
two other. To quantify this, we calculate the percentage of parti- 
cles lost between the 8th orbit (when the dust phase is incorpo- 
rated) and 104th orbit (the end of the simulation). Particles can 
be lost either by accretion onto the central star, accretion onto 
the planet or escape to large radii. The results of these tests are 
given in Table|2] We find that about twice as many particles are 
lost when (a,/3)=(l,2) than in the other cases. Simulations with 
(a,y6)=(0.1,0) and (0.1,0.5) give similar results, with the latter 
being a little more dissipative. The evolution of the disk with 
different SPH viscosity parameters is shown in Fig.[T] We can 
see that the (a,/3)=(l,2) simulations gives qualitatively differ- 
ent results, and that the (a,y6)=(0. 1,0) and (0. 1,0.5) simulations 
differ mainly in the corotation region and close to the edges of 
the gap where shock capture is more important. 

In order to keep the viscosity as low as possible, and con- 
sidering that results with /J = are not so different from those 
with p - 0.5, we keep a = 0.1 and /3 - for the remainder of 
our simulations. 



4. Simulation results 

Here we present the results from our suite of simulations, start- 
ing with our standard disk model, followed by our grain size 
tests and finally planetary mass tests. 



(11) 4.1 . Standard model 



In the standard model, we set up a 2.9 10 ^ Mq disk containing 
1% dust by mass made up of boulders 1 m in size. The disk 
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extends from 2 AU to 20 AU around a 1 Mq star and hosts a 
IMjup planet in a circular orbit at 5.2 AU. 

In Fig.[2]we show the top-down (x, y) and side-on (r, z) view 
of the density structure in the gas and dust phases of the stan- 
dard model. While the planet opens a shallow gap in the gas 
disk, the gap it creates in the dust layer is much more striking. 




r (AU) 



r (AU) 



Fig. 2. End state of the standard model after 104 orbits 
(1,230 years). Top: top-down view. Bottom; side-on view. Left: 
gas disk coloured by density. Right: dust disk coloured by den- 
sity. 

In Fig. |3]we plot the radial profile of the azimuthally aver- 
aged surface density of the standard model at five evolutionary 
stages (8.4P, 16P,32P,64P and 104P, where P is one period 
of the planet). Note that the dip in the radial density profile 
due to the planet induced gap is surrounded by two density en- 
hancements. The density enhancements are much sharper in 
the dust than in the gas even though the external enhancement 
gets broader over time due to feeding by external disk material. 
Particles that migrate inwards from the outer regions of the disk 
become trapped at the outer gap edge and over time more and 
more dust moves inwards and piles up at the outer gap edge, 
forcing this region of enhanced dust surface density to broaden. 
In the outer part of the disk, particles seem to be trapped by 



the 3: 2 resonance as noted bv lPaardekooper & Mellemal ( 12004 , 
I2OO6 '). We will return to this issue in Sect. 15.21 In the inner parts 
of the disk, on the other hand, particles flow through the reso- 
nances without becoming trapped due to the combined effects 
of the spiral waves created by the planet, the gas flow accreted 
by the central object, and the drag force that causes the dust 
to drift inwards faster than the gas once it has collapsed to the 
midplane. The internal over-density seen in Fig.[3]cannot be ex- 
plained by an accumulation of inwardly migrating grains, since 
they cannot cross the gap. It is instead due to the accumulation 
of dust at the gas pressure maximum. 

The gap is not exactly centred on the planet's orbit but it is 
shifted towards the inn er disk. This phenomenon was al so seen 
in the simulations of Paardekooper & Mellemal (2004) and is 
due to the fact that the dust does not directly feel the shock 
related to the creation of the gap. The shock only influences the 



dust via gas drag in the post-shock region. The flux of grains 
through the spiral wave is larger than if they were perfectly 
coupled to the gas and because the wave gradually depletes the 
gap, the density dip is deeper in the dust disk than in the gas. 

Figure H] shows the volume density of gas and dust as a 
function of distance to the star at the end of the simulation. The 
general radial dependence of the volume density is quite simi- 
lar to that of the surface density seen in Fig.|3] with similarities 
seen in the shape of the dust profile around the gap and the 
accumulation of dust at the corotation radius. Indeed, dust par- 
ticles tend to pile up where the gas concentrates, and Fig.|4]re- 
veals a denser region of gas at the position of the planet which 
is harder to see in the vertically integrated surface density in 
Fig. |3] but which was hinted at by the reduced scale height of 
the gas disk in the lower left panel of Fig.|2] The most striking 
feature of Fig.|4]is that, because of the very efficient vertical set- 
tling, the dust volume density reaches that of the gas in several 
regions of the disk, and even exceeds it at the inner edge of the 
gap where the dust piles up. Dust may therefore significantly 
affect the structure and evolution of the disk in that regime. 



g -13 
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Fig. 4. Volume density profile of the gas and dust of the stan- 
dard model after 104 orbits. The black points represent the gas 
particles and the coloured points represent the dust particles 
(coloured by density). Note that in this case the gas density has 
not been rescaled. 

In order to better understand the resulting structures in the 
dust disk, we ran a high resolution (A^~400,000) simulation of 
our standard model and studied the particle loss rate. The sim- 
ulations started with 200,000 gas particles and after 8 orbits 
the 200,000 dust particles were added. The particle numbers 
were compared at the end of the simulation. (Note that some 
gas particles are lost in the first eight orbits and hence we only 
compare particle numbers from the eighth orbit onwards.) In 
Table|3]we tally the number of gas and dust particles lost when 
R > Resc, when R - R^n (particles accreted by the planet), and 
when R < Rmnei (particles accreted by the star). 

We find that about 12% of the gas particles are lost while 
less than 1 % of the dust particles are lost. This is because dust 
particles tend to move inwards at the outer disk edge while gas 
particles tend to move outwards into the pressure void. At the 
edge of the planet gap, dust particles are trapped in the gas pres- 
sure maxima and so very few particles are actually captured by 
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Fig. 3. Azimuthally averaged surface density profiles for a gap created by a 1 Mjup planet at 5.2 AU in the standard model. The 
left panel shows the evolution of the 1 m sized grain dust layer after 8.4, 16, 32 and 64 and 104 planet orbits. The right panel 
shows the comparison of the gas and dust surface densities after 8.4 and 104 orbits. The gas density is scaled by 0.01 for direct 
comparison to the dust density. The vertical lines are the resonances as described in Fig.[T] 



Table 3. Particle loss rates from the 8th to the 104th orbit for 
high resolution runs of the standard MMSN model. Column 
one shows the accretion zone where particles are lost and 
columns two and three are the number of gas and dust parti- 
cles lost respectively (with the percentage in brackets). Row 
one shows the total number of particles (gas+dust) remaining 
at P = 104. 



Accretion 


A'gas A'dus, 


(Remaining) 


329,841 (87.5%) 


R > i?esc 


15,663 (4.2%) 202 (0.05%) 


R = ^crit 


18,746 (5%) 3,180(0.8%) 


R ^ dinner 


9,424(2.5%) 116(0.03%) 



the planet. At the inner disk edge, the gas particles again feel 
the pressure void and move inwards, while the dust particles 
remain where the gas pressure is highest. 

Figure |5] compares the resulting gas and dust surface den- 
sity profiles of the standard model for the low and high reso- 
lution simulations. The low resolution simulation captures the 
general shape and depth of the gap quite well, though it slightly 
overestimates the peak density at the gap edges and underesti- 
mates the dust resonant trapping by the planet. We are confi- 
dent that our low resolution simulations effectively capture the 
features of the planetary gap. 

4.2. Effect of grain size 

As discussed in Sect. [21 the scale height of the disk determines 
whether a planet of a given mass will open a gap in the disk, 
and from the results of BF05 we know that the evolution of 
the dust disk's scale height depends on the grain size. BF05 
found that large bodies (i > 1 m) are weakly coupled to the 
gas and essentially follow Keplerian orbits, while small bodies 
(s < lOyum) are very strongly coupled to the gas. In both cases, 
the dust disk remains flared but for different reasons. Only in 



— dust,hr 
— gas,hr 
dust,lr 
gas,lr 




Fig. 5. The resulting gap for the standard model for the low 
resolution (k: grey) and high resolution (hr: black) simulations. 
Dust profiles are shown as solid lines while gas profiles (scaled 
by 0.01 to compare with the dust) are shown as dashed lines. 



the inner part of the disk, where the gas density is highest and 
hence the gas drag is most efficient, does the disk flatten in 
the case of large grains. For intermediate sized grains (100 /vm 
< i < 10 cm) BF05 found that the drag strongly affects the dust 
dynamics, resulting in rapid settling to the midplane and sub- 
sequent strong inward radial migration. They found that 10 cm 
grains experienced the fastest radial migration. While the neb- 
ula conditions of BF05 are not identical to those used in this 
study, we expect the results to be qualitatively similar 

We ran a series of simulations to determine the effect of 
grain size on the induced planetary gap. Three grain sizes are 
tested: 1 cm, 10 cm and 1 m. Figs.|6]and|2]show how the shape 
and depth of the gap in the dust phase of the MMSN disk varies 
according to the dust grains size. (Note that the gas disk is un- 
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Fig. 6. Gap created in the dust disk for a range of grain sizes by a 1 Mjup planet at 5.2 AU after 104 orbits. Top panel shows the 
top-down view of disk and bottom panel shows the side-on view. From left to right: 1 cm, 10 cm and 1 m sized grains. The dust 
is coloured by density. 
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Fig. 7. Azimuthally averaged surface density profile for a gap 
created in the dust layer for different grain sizes (1 cm, 10 cm, 
1 m) by a 1 Mjup planet at 5.2 AU after 104 orbits. Gas sur- 
face density is scaled by 0.01 to compare with the dust and the 
vertical lines show the resonances as in Fig. [T] 



affected by a change in grain size and looks identical to the left 
hand panel of Fig. |2]) 

Similar to BF05, we find that the radial migration is most 
efficient for the 10 cm sized grains over most of the disk, which 
differs from the findings of Weidenschilling (1977). This dis- 
crepancy can be explained by the diff erent nebula parameters 
used, in particular the density. Indeed, Weidenschillinj ( 1977 ) 



noted that the radial migration velocity reaches a maximum for 
grains which have fs/forb - l/27r, where forb - 27r/QKep is the 
orbital period. Such grains have therefore Ts - I and, follow- 
ing Eq. (6), one can estimate the "optimal" grain size for radial 
migration at a given location in the disk by 



^opt 



Pdf^Kep 



(13) 



Figure |8] displays the variation of iopt with distance from the 
central star: decimetric grains are the most efficient radial mi- 
grators in the planet gap region, meter-sized grains migrate 
most efficiently from the outer gap edge out and decimetric 
grains take over again in the outer disk. Centimetre grains are 
only efficient at migrating in narrow regions on the very edges 
of the disk. 

As a result, while the 1 m boulders have the most rapid 
settling rate, the 10 cm grains have the most rapid radial mi- 
gration rate in the outer disk and produce a small, truncated, 
high density dust ring exterior to the planet. The density in- 
crease for the 10 cm grains in the external disk is compara- 
ble to that in the case of 1 m boulders due to two compet- 
ing effects: the settling is a little less efficient and leads to a 
thicker dust layer than in the 10 cm dust simulations, but more 
dust coming from the outer edge, where the radial migration is 
more efficient, accumulates at the external edge of the gap. The 
1 cm grains are most strongly coupled to the gas and thus their 
timescale of settling is much slower. This results in a broader, 
lower density dust disk. While the inner region of the 1 cm- 
sized grain disk appears evacuated of dust (Fig.|6]l, this is due 
to the smaller density increase around the gap (Fig. [T]) on the 
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Fig. 8. Size of the grains with the largest radial migration ve- 
locity at a given radius in the disk. 

timescales shown. This is further supported by the lower accre- 
tion rate seen in the 1 cm grains than in the 10 cm grains. 

As can be seen from Fig. |7] the gap gets deeper and wider 
as the gr ain size increases. Our find i ngs ar e in contradiction to 
those of IPaardekooper & Mellema ( 20061) who claim that the 
grain size has no influence on the width of the gap (though they 
do note that the gap gets deeper with increasing grain size as 
we do). While it is true that SPH will smooth the edges of the 
gap over 2h, we do not believe that this alone is the cause of the 
grain size dependent gap width. The h values are about 0. 1 AU 
around the gap edges in the 1 m and 10 cm cases and 0.2 AU 
in the noisier 1 cm case, which is far smaller than the smooth- 
ing of the gap edges seen in Fig. [7] We assume, therefore, that 
this discrepancy is due to the evolution of the dust layer's scale 
height with grain size which we can follow in our 3D simula- 
tions. Within this size range, the dust disk scale height dimin- 
ishes as the grain size increases, and hence the planet opens a 
gap more rapidly in larger grain sized dust layers. 

For the 1 m boulders case, particles are found in corotation 
with the planet. We would expect that some particles would be 
found at corotation for all grain sizes, but we assume that this is 
only seen in the 1 m case because the gap is deeper and wider 
than in the other two cases, making those particles caught in 
cor otation eas i er to s ee. 

Rice et al. I (l2006h find that the pressure gradient at the outer 



gap edge act as a filter, allowing smaller particles to penetrate 
into the inner disk while holding back larger grains. Their tran- 
sition size is 5 < lOfim or smaller, and depends little on disk 
parameters. The dust grains we use in our simulations are all 
larger than this critical size and we indeed find that particles 
outside the gap do not migrate inwards, leaving the inner disk 
to drain by accretion on the central star 

4.3. Effect of the planetary mass 

We also investigate the effect of the planetary mass on the evo- 
lution of the gap in the dust disk. We return to our standard 



MMSN model which contains 1 m sized grains and vary Mp 
from 0.05 Mjup to 1 Mjup. The resulting disk morphologies are 
shown in Fig. |9l (Note that in this series of tests, the increas- 
ing planet mass does affect the gas disk distribution as well. 
However, the effect is quite small and so we do not include gas 
distribution figures.) 

The azimuthally averaged surface density profiles after 104 
orbits are shown in Fig.[TO]for the dust and gas phases. An 0.05 
Mjup planet produces a decrease in the dust surface density in 
the vicinity of its orbit, whereas the gas surface density is not 
affected until the planet is more massive than 0.2 Mjup. 

As expected, the gap gets deeper and wider as the mass of 
the planet increases. In the inner disk, the dust surface den- 
sity increase is due to the dust accumulation at the gas pressure 
maximum. While grains cannot cross the gap from the outer 
disk, in the inner disk grains move towards the pressure maxi- 
mum coming from both sides of the maxima. The peak of the 
dust surface density is dictated by the efficiency of this process 
and does not depend on the gap depth. However the density 
peak is broader for smaller mass planets because of the coiTe- 
spondingly shallower gap. For low mass planets, material can 
still flow from the outer to inner disk and therefore we have 
more mass in the inner disk than when the gap is deep enough 
to prevent this inflow. In the case of larger mass planets, the 
density peak gets thinner before its maximal value starts to de- 
crease. 

In the outer disk, the maximum surface density increases a 
little with the mass of the planet but reaches almost the same 
value for the 0.2, 0.5 and 1 Mjup planets. The peak is almost 
as broad in each case but is shifted to larger radii as the gap 
widens with increasing planetary mass. 

For the two most massive planets modelled (Mp = 0.5 and 1 
Mjup) some dust particles appear to become trapped near coro- 
tation. Once again these are the simulations in which the gap 
is deepest and widest and hence it is easier to see particles in 
corotation. 

In all cases but for the 1 Mjup planet, the gap is again shifted 
sli ghtly towards smaller radii t han the planet orbit. 

Paardekooper & Mellemal ( 20061) find that the smallest 



planet mass able to open a gap in a dust disk comprising 
1 mm-sized grains located at 13 AU is 0.05Mjup for simu- 
lations of a MMSN disk. While their nebula parameters are 
not identical to ours and their models are 2D and less viscous 
(ffss - 0.004), there is still general agreement between our re- 
sults. From Eq. ([T]i we would expect a larger mass planet is 
required to open a gap at smaller radii, but with effective dust 
settling, the disk thickness is small at small radii which reduces 
the planet mass required to open a gap. Thus we also find that 
an 0.05Mjup planet at 5 AU (for our disk parameters) can open 
a gap in the dust disk. 



5. Discussion 

Structures created by planets in dusty disks are more diverse 
than those created in the gaseous disks. With only aerodynamic 
drag, we already have the possibi lity of creating either a disk 
with a large central hole or a ring. iRice et alj ( l2006t) also found 
that, for certain grain sizes, planets can create a central hole in 
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Fig. 9. Gap created in the dust disk comprising 1 m grains by a planet of varying mass at 5.2 AU after 104 orbits. Top panel 
shows the top-down view of disk and bottom panel shows the side-on view. From left to right: a planet of mass 0.05, 0.1, 0.2, 0.5 
and 1.0 Mjup. The dust is coloured by density. 
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Fig. 10. Azimuthally averaged surface density profile for a gap created in the disk after 104 orbits by planets of various masses 
after 104 orbits. The left panel shows the evolution of the dust layer containing 1 m sized grains and the right panel shown 
the evolution of the gas layer which has been scaled by 0.01 for direct comparison with the dust. The vertical lines show the 
resonances as in Fig.[T] 



the disk. Other processes can f orm similar structures, e . g. disk 
clearing by photoevaporation (' Alexander & Armitage 2007 ). 
The influence of radiation pressure and Poynting-Robertson 
drag would clearly lead to other more detailed structu res for 
smaller grains (see, e.g. Takeuchi & ArtymowiczlEoOl ). 



5.1. The effect of the inner boundary 

Recently ICrida & MorbideUil (l2007h showed that the location 



of the inner boundary in planetary gap simulations can have a 
dramatic effect on the resulting structure of the innermost part 
of the disk. They found that a planet will open an inner cavity 
(i.e. an inner hole) in the disk if the planet is more massive that 
the disk or if the inner disk edge, rin, is too close to the planet. 
They show that the disk surface density profile in the inner disk 
is lowered by a factor 1 - -<Jrijr^ in the presence of a planet. 



We re-ran our standard model with the inner boundary of 
the disk at ri„ =0.1 and 0.25 to compare the resulting inner 
di sk structure with our Q,, j^O .4 simulations. Using the formula 
of lCrida & MorbideUil (l2007h . we should find that the inner sur- 
face density drops by 32% when rin = 0.1, compared to a 63% 
decrease when rin - 0.4, wh ich is what is seen in Fig. fm (and 
compares well with Fig. 7 of lCrida & Morbidelhll2007h . 

Crida & MorbideUil tOOl h speculate about what would 



happen to dust in this inner disk region. They suggest that the 
dust would accumulate at the relative maximum of the gas den- 
sity, which is indeed what we find, and that a thin dust ring 
would result, just as we see in our disks in Fig.|6| 

While it is clear that the location of the inner boundary has a 
dramatic effect on the density distribution of the inner disk, the 
results we have presented for different grain sizes and different 
planet masses should not change for a fixed ri„. 
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Fig. 11. Surface density profile of the gas in the standard model 
for three different inner boundaries; rin = 2.08 AU, 1 .3 AU and 
0.52 AU. 




Fig. 12. Eccentricity (e) versus semi-major axis (a) for the stan- 
dard MMSN disk with a 1 Mjup planet orbiting at 5.2 AU after 
104 orbits. Left panel: without gas drag. Right panel: gas drag 
included. 

5.2. Trapping in the 3:2 resonance 

Paardekooper & Mellemal (l2006l) suggested that particles be- 
come trapped in the 3:2 external resonance with the planet. 
While we see a clear density increase in the vicinity of the 
3:2 resonance of our standard MMSN disk in Fig.|3] it appears 
to be too broad to caused by resonant trapping. To investigate 
further, we plot the eccentricity of the dust grains versus semi- 
major axis in Fig. [12] In these kind of diagrams, resonances 
appear as thin vertical lines and as a V-shaped pattern at the 
edges of the gap. We ran a simulation without the drag force 
in which dust particles are fully decoupled from the gas. The 
results are shown on the left panel of Fig.[T2]and the signatures 
of resonant trapping can be clearly seen. The rather high eccen- 
tricity of the dust particles is due to the initial conditions - the 
dust was initially given the velocity field of the gas and without 
drag there is no mechanism to damp the gas velocity disper- 
sion. In the right panel of Fig. [121 we see that when simulations 
are run with gas drag, the drag is very efficient at damping the 
particles' eccentricity and consequently all resonant signatures 
disappear. Furthermore, if we compare our results for different 



grain sizes (see Fig.|2ll and for different planet masses (Fig.llOll. 
we see that the particle pile up at the outer edge of the gap does 
not always coincide with the 3:2 resonance for the same evo- 
lutionary phase (104 planet orbits). This seems to depend quite 
sensitively on the exact disk parameters. 

We conclude, therefore, that the accu mulation o f dust 
that we and other authors ( Paardekooper & Melleni^ 12006 ; 
Alexander & Armitagel2007 ) notice close to the outer gap edge 



is not due to resonant trapping. Resonant trapping will depend 
on the efficiency of the gas drag for the grain sizes under con- 
sideration. 



5.3. Observations of protoplanetary d/sks 

Our results have implications for recent predictions of obser- 
vations of protop lanetary disks host i ng pla net gaps. As men 



tioned in Sect. 1. IWolf & D'Angelol (l2005h and lVarniere et al 



(l2006h lave used the results of 2D hydrodynamics simulations 
to produce synthe tic images of protoplanetary disks following 
the suggestion of [Paardekooper & Mellemal (l2004l) that such 
disks would be visible with ALMA. Such 2D simulations as- 
sume that the gas and dust are well mixed, which our results 
clearly demonstrate is not the case. 

The results presented here cannot directly be used to con- 
struct synthetic observations because we only consider one dis- 
crete grain size per simulation. However, since we don't con- 
sider interactions between dust grains and since the structure of 
the gas disk doesn't vary with grain size, we can stack our re- 
sults for different grain sizes to reconstruct a grain size distribu- 
tion. With a finer size sampling than presented here, we will be 



able to create synthetic scattered light images (see iPinte et al 



2007h . 

Because our general findings show than the gap is generally 
more striking in the dust disk, we suggest that predictions of 
observations of protoplanetary disks are too pessimistic. Our 
results show that the density contrast around the gap can be 
very strong - up to four orders of magnitude (see Fig. [Hi - and 
this would be detectable. Our results, together with those of 
IPaardekooper & Mellemal ( l2006l) . show that even if one varies 
the nebula parameters, the result is robust: an 0.05Mjud pl anet 
o pens a gap in the dust disk. As suggested by Varniere e t al.l 
(120061) the resulting strong density enhancements at the edges 
of the gap should be easily visible with ALMA. 

While the interactions of the dust with the gas disk and 
the planet are complex, our simulations can help constrain the 
planet mass and grain size distribution in future high angular 
resolution observations of protoplanetary disks. 

5.4. Enhanced planetesimal growth 

We have shown that the presence of the planet triggers an accu- 
mulation of grains at the external edge of the planetary gap, and 
the resulting density enhancement may well favour the growth 
of planetesimals in the region. This can be put in perspective 
with the results of Tsig anis et al. j2005 '), who propose a sce- 
nario for the formation of the Solar System that reproduces the 
characteristics of the giant planets as well as that of the Trojans 
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and explain the late heavy bombardment of the Moon. This 
scenario relies on the controversial hypothesis that Jupiter and 
Saturn formed closer to one another than they are today and 
subsequently migrated, one inwards and the other outwards. 
While our results are rather extreme in that one of the planets 
is fully formed whereas the other is only at the planetesimal 
stage, even a planet embryo (well before it could carve out a 
gap) would create regions of pressure enhancement in its vicin- 
ity where dust would accumulate. Such dense rings would be 
ideal sites for the concentration of solid pa rticle s, and may lead 
to faster growth of plane tesimals (jPurisen et al.ll2007l) . 

5.5. Planet migration 

In this paper, we do not investigate the influence of planetary 
migration. The planet remains on a fixed circular orbit in our 
simulations and momentum exchange between planet and disk 
is not taken into account. We would expect, however, planet mi- 
gration to have an effect on the dusty gap formation. Because 
the total dust mass is two orders of magnitude lower than the 
gas, one would expect the migration to be driven by gas and not 
dust. However, would this really be the case? Ignoring the in- 
fluence of turbulence, once the dust has settled to the midplane 
its density tends to be comparable to that of the gas (see Fig.|4|i. 
One may speculate that the planet might be locked to the dust 
evolution and gas drag imphes that dust flows faster towards 
the central object than gas. However, turbulence would need to 
be taken into account to fully study such a scenario. We will 
investigate the influence of migration in a dusty protoplanetary 
disk in a future paper 

6. Conclusions 

We have conducted a series of 3D numerical simulations of 
two-phase (dust-Hgas) protoplanetary disks to study the be- 
haviour of the dust in the presence of a planet. We ran a se- 
ries of experiments with a minimum mass solar nebula disk, 
varying the grain size and the planet mass. 

We find that gap formation is more rapid and striking in the 
dust layer than in the gas layer Varying the grain size alone 
results in a variety of different structures which may well be 
detectable with ALMA. For low mass planets in our MMSN 
disk, a gap was found to open in the dust layer while not in the 
gas layer 

Grains accumulate at the external edge of the gap but we 
find that this is not due to resonant trapping when gas drag is in- 
cluded. This is more likely due to the dust concentrating at the 
gas pressure maxima. Contrary to previous 2D work, our 3D 
simulations demonstrate that there is a dependence of the width 
and depth of the gap on grain size. Simulations like these can 
be used to help interpret observations to constrain the planet 
mass and grain sizes in protoplanetary disks. 
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